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Abstract. The kinetics of the long-range spherical model evolving from various 
initial states is studied. In particular, the large-time auto-correlation and auto- 
response functions are obtained, for classes of long-range correlated initial states, and 
for magnetized initial states. The ageing exponents can depend on certain qualitative 
features of initial states. We explicitly find the conditions for the system to cross over 
from ageing classes that depend on initial conditions to those that do not. 



1. Introduction 

Slowly-relaxing thermodynamic systems are being intensely studied for decades now, 
and for numerous reasons. One of the main objectives is to understand genuine 
non-equilibrium states: effecting their description, finding and categorizing the large- 
scale properties, identifying the necessary qualitative features of each universality class, 
determining sufficient exponents that would enable us to label the classes, and answering 
a host of other similar questions usually posed in equilibrium statistical physics. 

Systems that exhibit phase transitions relax slowly and may not equilibriate when 
either quenched to critical temperature or below. For critical quench the relaxation 
is inhibited due to long-range correlations [I], while when quenched from above to 
below the critical temperature, it is due to the competition between locally-equilibriated 
domains for a global equilibrium state [2]. In either of these cases, after waiting for a 
long time t w , much longer than the transient microscopic time-scales, these systems 
which are not in equilibrium are expected to become statistically scale invariant, when 
scaled by a characteristic time-dependent length-scale L(t). The length grows in time 
algebrically, L(t) ~ t 1 ^, defining the dynamic exponent z. The scaling hypothesis, 
if applicable, proposes that the two-point functions of observables, say correlation- or 
response- function, take the form 



for large observation time t > t w ; where Cj(x, t) refers to some local observable at time 
t and position universal exponent, i*Yj is a scaling function. Here it is assumed 
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that the system is translational- and rotational- invariant. Furthermore, in the ageing 
regime, namely when t 3> t w , the asymptotic behavior of the function is as follows, 



and defines an ageing exponent A^-. The exponent though may not be an independent 
exponent. For instance, when t Pd t w , and |x — y| <C L(t), and further under the 
conditions that would establish equilibrium correlations of size L(t) in time t, the two- 
point function of observables at critical quench would behave as (Cj(x, t)Oj(y, t w )) ~ 
|x — y| _Aij , where Ajj is an equilibrium exponent. If this is so for (nearly) equal-time 
observables, then consistency with scaling hypothesis would imply = Aij/z, and 
hence is not an independent exponent. 

Scaling hypothesis holds in many slowly-relaxing systems, and is tested 
experimentally and theoretically [3.]. Long-range spherical model is one of the few models 
that can be solved exactly. This model when evolved from an uncorrelated initial state 
shows scaling behaviour at large times Initial conditions are known to affect the 
large-scale properties in some short-range models[5], El [TJ [HJ El EH1 EH E2l H3J- So it 
would be interesting to know how, in general, the initial conditions influence the scaling 
behaviour, and in particular, that of long-range models. It is this issue that we address 
in this paper, after having solved the non-conserved long-range spherical model evolving 
from different initial states. We evaluate the correlation and the response functions for 
various initial conditions, and determine their scaling behaviour, and also find how the 
ageing exponents depend on certain features of the initial state. In some cases the 
aij exponents turn out to be independent, implying that the spacial correlations within 
regions of size L(t) after a time t are not governed by canonical equilibrium fluctuations. 

The layout of the paper is as follows. In the next section, the model is defined 
and an outline of the formal solution is provided. In section [3X the scaling behaviour 
of auto-correlation and auto-response functions, with long-range initial correlations is 
investigated. In section HJ the model is analysed with an initial state that has no 
correlations but has non-vanishing magnetization. Section [5] concludes with a brief 
summary and a few remarks. 

2. Long-range spherical model 

A spherical model is exactly solvable not only in equilibrium but also when evolving 
with relaxational dynamics [14]. The techniques used in solving that model can be easily 
extended to a long-range spherical modeljl]. In this section, we first define the long- 
range model and then solve for the correlation and response functions. 

2.1. The model 

The equilibrium long-range spherical model [15] is described by the Hamiltonian 




(2) 




(3) 
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where the coupling between the spins is long range and given by 

^,o = ^ixr (d+CT) , (4) 

such that J2x ^x,o = 1, which in turn fixes the constant J c . The variable S x respects the 
spherical constraint Y[S] = 0, where 

nsi = ^ j>(*) 2 - 1, (5) 

xeA 

and, x runs over the sites of the ci-dimensional lattice A, and iV is the total number of 
sites on the lattice. 

The coarse-grain dynamics near the critical temperature or the coarsening dynamics 
when quenched from above to below T c , for the non-conserved order-parameter S x is 
given by following Langevin equation, 

5H 
~~5S X 

(Vx(t)Vx>(t')) = 2T8(t-t')5 x , xl , (7) 

where Z is the Langrange multiplier that is determined by the constraint. Assigning 
Stratonovich convention to the stochastic dynamics will lead the constraint to the 
expression, 



dtS x {t) 



-S x {t)Z(t)+r) x (t), (6) 



(8) 



2.2. Formal solution 



We will first separate the Gaussian and non-Gaussian terms in equation ([6]) so as to 
solve. We will see that the fluctuations of the Lagrange multiplier can be neglected 
when restricted to local observables in the limit A" — > oo. Furthermore, the spherical 
constraint can be replaced by a mean-spherical constraint in this thermodynamic limit. 
We will then explicitly write down the solution in the Gaussian regime using this mean- 
spherical constraint. 

2.2.1. Fluctuations of the Lagrange multiplier 

The fluctuating part of Z(t) is of O (l/VN^j, and hence S x can be written as, 

S x (t) = Sx (t) + J=f x (t), (9) 



where s x captures only the Gaussian fluctuations while £ x accounts for the remaining 
non-Gaussian fluctuations. 

Within the Gaussian approximation the mean-square fluctuations in Y[s] can be 
estimated from the following relation, obtained using Wick's theorem and translation 
invariance: 

(? 2 m>-(^]} 2 4e^(m), (io) 
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where C X) o(£, i) is the equal-time spin-spin correlation function. In the case of 
uncorrelated initial state, the late-time behaviour of the correlation function is as given 
below [I], 

r (tt\ J G (M* -1 '*) > T <T C 
Cx ' o(M) \ 1*-*/°G (\x\t-V') , T = T C ' 

where, G(x) oc J k exp (ik • x — 2|k| CT ). Hence we obtain, 

(?•[.]>-<?[.]>• «{^_ (12) 

Therefore, in the limit N ^ oo these fluctuations become negligible for late-times 
t < t* ~ N a / d [16], both at and below the critical temperature T c . Even at T c the 
threshold time, t*, is unchanged since the phase transition occurs only for < o < d. 
Thus, taking first the N —>■ oo limit is a necessary condition to have the same asymptotic 
behaviour for both spherical and mean-spherical models in the Gaussian regime. 

Substituting the decomposition of spin variables, as given in equation (jHJ), into 
equation (JSJ) and keeping only the leading subdominant terms yields the expression 



Y[S) « (Y[s)} + i 



(13) 



Here the non-Gaussian term is of the same order in iV as the Gaussian fluctuations, and 
hence can play a role in certain global observables. 

2.2.2. Solution in the Gaussian regime 

The Gaussian approximation of the Langevin equation ([6]) in the Fourier space is given 
as 

d t s k {t) = -(u h + Z(t)}s k (t)+Ti k {t), (14) 

where u k is the Fourier transform of —J x ^, and the noise T] k (t) has the variance, 

(Vk{t)Vk'{t')) = 2T5(t - t')N5 k+k ,. (15) 

In the continuum limit, N5 k gets replaced by (2n) d 5(k), and iV -1 J^ fc is replaced by 
(2ir)~ d J d d k. In this limit, for small k, the function uo k — > Cq + ci\k\ a + C2\k\ 2 + 0(|/c| 4 ), 
where c n are lattice-dependent constants. We shall restrict a < 2, for a > 2 the relevant 
model is the usual short-range spherical model. Further, without any loss of generality 
we shall set 5q = by absorbing it into Z(t). 
The solution of the above equation is 



s k {t) 



s k (0)+ / dTe^ T ^/g^) Vk (r) 



(16) 



where the yet to be determined constraint function g(t) = exp(2 J Q * drZ(r)). It is often 
convenient to write the above equation as 

s k {t) = R k {t,0)s k {0)+ [ drR k (t,r)r] h (r), (17) 
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where Green's function R k (t,t') satisfies the equation, 

(d t + u k + Z(t)^R k (t,t') = 6{t-t'), (18) 

with the initial condition R k (t,t') = for t < t! . It is easy to see that R k (t,t') is also 
the Fourier transform of the response function R x fl(t,t'), and is explicitly given as 



R k (t,t') = e-"^J^6(t-t'), (19) 

where 9(t) is the step function. 

If the initial conditions for the spins s x are chosen to be translation invariant with 
a variance, 

MOMO)) = N6 k+k ,C k (0,0). (20) 
then, using equation (fT6l) . we obtain the correlation function 

(s k (t)s k ,{t')) = N6 k+k ,C k (t,t'), (21) 

where 

C k (t,t') = 7 gM e-^- f, 'C fc (t m ,t m ), (22) 
and t m = min(t,t'), and the equal-time correlation function 



C k (t,t) 



C fc (0, 0) + 2T / dre 2 ^ T g(r 



t 



(23) 



git) 

Note that C k (t,t) is also the Fourier transform of the correlation function C x (t,t) = 
(S x+y (t)S x (t)) . We now impose the mean-spherical constraint, 0^[s]) = 0, which then 
implies that the following condition should hold at all times, 

I^C fc (M) = l- (24) 

k 

This condition not only fixes g(t) but also restricts the choice of initial conditions. 
Substituting expression (T23]) in the above equation gives g{t) to be the solution of the 
following Volterra equation, 

g(t) = A{t) + 2T f drf(t - r)g(r), (25) 



with g(0) = 1, where 

Mt) = ^J2 e ^ ktC ^ ^ ( 26 ) 



N 
k 



and 



N 

k 



Note that A(0) = 1, due to spherical constraint. The function A(t) depends on the 
initial correlations, and hence g(t). We shall solve equation fT25l) and determine g(t) for 
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various initial conditions in latter sections, and explicitly evaluate the correlation and 
response functions. 

The non-Gaussian variable £ x (^) can be completely expressed in terms of Gaussian 
observables as has been done recently for the short-range spherical model [10]. But 
we shall not express this here, since we are interested in the ageing properties of local 
observables, for which the non-Gaussian corrections vanish in the limit N — > oo. 



3. Long-range correlated initial conditions 

In this section, we will evaluate the auto-correlation and auto-response functions, when 
evolved from a long-range-correlated initial state. To this end, we will first calculate the 
constraint function. 



3. 1 . Constraint function 

The large-time solution of the equation (125]) for the constraint function g(t) not only 
depends on the temperature T but also on the form of the initial correlation function 
Cfc(0, 0). The case with uncorrelated initial conditions in this model was studied 
recently [4]. We now investigate the behaviour when evolved from a class of initial 
states with long-range spin-spin correlations. In other words, we choose (7^.(0,0) such 
that its small k behaviour is given by 

C fc (0,0)^c \k\°°, (28) 



and also satisfies equation ( 1241) . In the real space, the large- distance initial correlations 
become C x (0, 0) ~ c' \x\~( d+ao \ The exponent o~o is restricted by the condition 
d+o-Q > 0, on the grounds that the correlations should decrease with increase in distance. 
Short-range spherical model with similar initial state was also studied recently [8J. 
The Laplace transform of equation (125|) results in the expression, 

»W = i-2T /l(P )' (29) 

where, gi(p), Al(p) and flip) are the Laplace transforms of g(t), A(t) and f(t), 
respectively. The small p behaviour of gh{p) is sufficient for deducing the large-t 
properties of g(t). The function flip), for small p, is given as 

oo 

hip) = -A p- 1+d /° + A n+1 (-p) n (30) 

n=0 

where A = |T(1 - d/a)\A , and A n = N" 1 ^ k {2uj k )~ n - j k (2cx|A;| CT )-", and A n = 
J k exp(— 2ci|/c| CT )|co/c| n<J °. Similarly, Ai(p) is given by 

oo 

A L {p) = -B p- 1+ ^V° + B n+1 (-p) n (31) 

n=0 

where, B = \T(1 - (d + <t q ) / <t)\coA x , and B n = N' 1 J2 k (2co k )- n C k (0,0) - 



The role of initial conditions in the ageing of the long-range spherical model 



7 



Regime 


Conditions 


g (a, (T , d) 


ip(a, a , d) 


I 


d/2 < o < d, 
—d < cr < o — d 


AiB 
y—ao/a) 


-1 - oq/o 


II 


< a < d/2, 
—d < cr < o — d 


A 1 B Q 
A 2 l (2—(d+ao)/u) 


1 - (d+ (T )/(T 


III 


d/2 < o < d, 
a > a — d 


a t~v 1 i ^ / ~\ 
Aol (, — l+a/cr; 


-2 + d/cr 


IV 


< a < d/2, 
o"o > o — d, do > —a 


A 2 





V 


< (T < d/2, 

(Tq > (T — d, (Tq < —O 


A 1 B 1 
A 2 






Table 1. Values of go and V: as defined in equation (|32|) . characterizing the 
asymptotics of <7(i) in various critical regimes. 



There is a phase transition at temperature T c = 1/2 A\. At the critical temperature 
T c , the small p behaviour of giip) can easily be obtained by substituting equations (1301) 
and (13TT) into (1291) . In this limit, (?lQ») ~ fi'orX'*/' + l)p _ ^ +1 \ and hence the large-time 
form of g(t) is given by 

g(t)^g (a,a ,d)t^' d \ (32) 

where the constant go(o~, cr , d) and the exponent i/)(cr, o~o,d) depend on the values of 
a and cr , for a given d There are four regimes, as specified in the table([T]), each 
distinguished by the form of ip. Though the form of g(t) in case-IV and case-V is the 
same, they are listed separately for, as we shall see later, they get distinguished by 
ageing exponents and the fluctuation-dissipation ratio. 

For T < T c , in the small p region the function gi(p) ~ (1 — T/T^^Al^p), and 
hence 

g(t)*A 1 (l-pj r^/* (33) 
3.2. Auto- correlation and auto-response functions 

We first write the auto-correlation and auto-response functions in terms of g{t), A(t) 
and f(t), and then analyse their late-time behaviour. The correlation function is given 
by 

C x {t,t w ) := (S x+y (t)S y {t w )) = l^ e ifc - x C fc (t,U, (34) 

k 

and hence the expression for the auto-correlation function, C{t, t w ) := C x= o(t, t w ), using 
equations ff22l) and ff23l . becomes 



C(t, t h 



A ( t -^-)+2T J\rf { t A^-r)g{r 



.(35) 
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The spin response to magnetic field is obtained by perturbing the Hamiltonian, 
H(t) — > H(t) — f dT J £2 x S x (r)h K (T), with an impulse of magnetic field at time t w , and 
is given by 

*<W*)>* 

h=0 



5h y (t WJ 



= ^^'^f^ (36) 

for t > t w , and zero otherwise. The auto-response function, R(t,t w ) := R x =o{t,t w ), 
therefore takes the form 




R(t,t w ) = J^f — ^ . (37) 



3.2.1. Critical quench 

The asymptotic behaviour of f(t) ~ AQt~ d ^ a , and A{t) ~ Ait~( d+a °^ CT , is independent 
of temperature. Making use of these asymptotics and that of g(t) at T c as given in 
equation ([3"2l . the expression fl33|) at large t and t^, reduces to 

C(t, t w ) = 2 d ^A t 1 -^ [t w {1+ ^ /a) F lc (x) + F 2C (x)] , (38) 

where, x = t/t w and 

F lc (x) = -4^2< Tt,/,J x~^ 2 (x + i)~( d +°o)/^ (39) 

F 2C {x) = 2T C x^' 2 [ dv{x + 1 - 2vY dt(J v^. (40) 
Jo 

Thus the auto-correlation function has the scaling form 

C(t,t w ) = 2 d ^A tZ b f c (x), (41) 
with the following three possibilities: 

(i) b = -1 + d/a and f c (x) = F ic (x) + F 2C (x), when 1 + if) + a /a = 0; 

(ii) b = — 1 + d/a and /c^c) = F 2C (x), when l + ?/>-|-cro/cr>0; 

(iii) 6 = if) + (d+ do)/ " an d fc — Fic( x ), when 1 + if) + (T /cr < 0. 

One of these possibilities occurs in each of the five regimes mentioned in table([T]). In 
the limit x — >• oo 

fc(x) ~ oT A ^, (42) 

where the dynamic exponent 2 = 0". The values of 6 and Ac in various regimes are 
listed explicitly in table([2]), with one fine-tuned exception in regime-I. This exceptional 
case is when not only 00 = —0 but also the amplitude of the initial correlation is tuned 
such that it exactly corresponds to correlations of the equilibrium long-range spherical 
model. In this fine-tuned case, the system is always at equilibrium and the value of the 
ageing exponent Ac = d. 
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Similarly, using the asymptotics of f(t) and g(t) in equation (137]) . we obtain the 
expression for the response function for large times as 

R(t, t w ) = 2 d ^A t- d ^x-^ 2 (x - l)- d/a , (43) 

where x^l. Here again the auto-response function takes the expected scaling form 

R(t,t w ) = 2 d /°A t~ 1 - a f R (x), (44) 

where a = — 1 + d/a and fii(x) = x~^^ 2 (x — l)~ d /°", which in the limit x — > oo behaves 
as 

f R (x) ~ oT A ^, (45) 

with \r = d + ai^/ 2. The values of 6 and Ar in all the five regimes are listed explicitly 
in table©. 

Another interesting quantity to look at is the fluctuation-dissipation ratio (FDR). 
A deviation of this ratio from unity indicates that the state is nonequilibrium. FDR, 
which is defined as 

* (t < tJ:=r ror (46) 

can be obtained explicitly up on using equations (!38|) and (|43l) . and is given by the 
expression 



x 



-V>/2 



(x - l)- d/a 



X(t, t w ) = T c — — — — -A= '- = , (47) 

t -^^F lc (x) + F 2C (xy 

where, x ^ 1, and 

F ic {x) = -U+ F lc (x) - x^F lc (x), (48) 

\ a J dx 

F 2C (x) = (l - ~) F 2C (x) - x^F 2C (x). (49) 

The large x behaviour of F ic (x), using equations fl39|) and (jUJ, is easily obtainable and 
is given by 

F ic (x) ~ -^2 CT °/ CT (-k )x- {kl+id+ao)/a) , (50) 

where either k$ = ki = ip/2 when ?/> 7^ 0, or A; = (<i + o~ )/a and &i = 1 when ip = 0. 
Similarly, the large x behaviour of F 2 c(x), using equations (1401) and (]49l . for t/> 7^ is 
given by 

F 2C (x) ~ T c 2 -±tx-^ +d /°\ (51) 
1 + if) 

while for tp = 0, 

F 2C (x) = T c [(x + l)- d ^ + (x- !)-*'*] . (52) 
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Regime 


b 


Ac 


a 


Ar 




I 


-1 + d/a 


d+ (a -a)/2* 


-1 + d/a 


d- (a + a )/2 


0<r+cr ,0 


II 


1 


(d + a + a )/2 


-1+d/a 


(d + a-a )/2 





III 


-1 + d/a 


3d/2 - a 


-1 + d/a 


3d/2 - a 


1 -a/d 


IV 


-1 + d/a 


d 


-1 + d/a 


d 


1/2 


V 


(d + a )/a 


d + a 


-1 + d/a 


d 






Table 2. Ageing exponents and asymptotic FDR in various critical regimes; *except 
for the fine-tuned case, where Ac = d 



Using the above expressions for Fic{x) and F 2 c(x) in equation f HTj) . it is easy to evaluate 
the asymptotic FDR 

X°° := lim (]im X{t,t w )) = lim ( lim X{t,t w )L =t/s ] . (53) 

The asymptotic FDR in all cases is obtained as specified in table(T5]). 

A few remarks on the ageing exponents follow. Both regimes III and IV are 
independent of initial conditions, and the exponents are same as in the case of the 
uncorrelated initial state [I] . Note that in both these regimes the condition a + ao > 
is satisfied. A physical interpretation of this inequality is that the initial correlations 
are weaker than the equilibrium correlations towards which the system is evolving. This 
is easily deduced by comparing equilibrium correlations C eq (x) ~ |x|~( d ~°") [15] with 
initial correlations Cj n (x) ~ |x|~( d+cr °- ) . Also note that the auto-correlation exponent 
b in the case of uncorrelated initial state 0] is not independent, since we can obtain 
C eq (x) ~ |x| _bz , when b = — 1 + d/a is chosen. 

In regimes II and V, a + ao < 0, and hence can be expected that the ageing 
exponents will be influenced by the initial conditions, which indeed is the case. 

Regime I is where a + o"o can be positive, or negative, or even zero. Here some of 
the exponents, namely a and b, are blind to initial conditions while others, namely A<7 
and \r, are not. The earlier-mentioned physical picture seems to fail in this regime. 
When a + a = 0, the asymptotic FDR X°° = 1, though in general Ac ^ Xr- The 
only exception is the fine-tuned case where the system is at equilibrium to begin with. 
In all the regimes which depend on initial conditions, Ac < \r, suggesting that the 
correlations decay slower than the responses. 



3.2.2. Phase ordering 

We shall briefly discuss the scaling behaviour when the system evolves from the long- 
range initial conditions below critical temperature. Substituting equation fl33|) into 
expression (1351) . results in the following late-time auto-correlation function when T < T c : 

T\ {( x + lY\-( d+ ^l 2 ° 



^Hi-^j(V) * (54) 
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Thus scaling behaviour is exhibited even while coarsening; with the exponents 6 = 
and Xc — (d+ o~ )/2. 

Similarly, substituting equation ( 1331) into expression ( 1371) gives the late-time auto- 
response function as 

R{t, t w ) = 2 d/a A t- d/cT x {d+ao)l2a {x - l)- d/a , (55) 

and which shows scaling behaviour with the exponents a = —1+d/ a and Xr = (d—cr Q ) /2. 
It is easy to see, from equations (1541) and ( [551) . that in the phase-ordering regime the 
asymptotic FDR vanishes, or rather X(t,t w ) — > as t w — > oo for any fixed x. 

4. Magnetized initial conditions 

In this section, we evaluate the late-time auto-correlation and -response functions, within 
the Gaussian approximation, when the system evolves at critical temperature from an 
initial state having a non-zero magnetization, 

^5>*(°)> = ^(sk(0))\ k=0 = m , (56) 

X 

and uncorrelated spin-spin correlations, 

(s x (0)s y (0)) = a-ml)5 x , y + ml (57) 

which respect the spherical constraint. 

Recent studies in short-range spherical model [TOl [13] have shown that the non- 
Gaussian corrections are not relevant for local observables, and only become significant 
for global observables, when magnetized initial conditions are considered. These results 
will also hold for the long-range spherical model, with magnetized initial conditions. 

4-1. Constraint function 

Equation (1571) implies (7^(0,0) = (1 — ttIq) + m^NSkfi, which upon substituting into 
equation (1261) gives 

A(t) = (l-m 2 )f(t)+ml. (58) 

Using the Laplace transform of the above function in equation (|29|) . and then performing 
an inverse-Laplace transform in the small p region, gives the following late-time 
behaviour of the constraint function, 

g(t)*g M (a,<£)r™(l + ^-), (59) 

where tu = (1 — 20)(1 — ml)/2T c ml, is a time scale associated with the initial 
magnetization, and 

/ 0, 0<a<d/2 

II- d/2a, d/2<a<d ' 1 ' 

, I -ml \ 1/A 2 , 0<a<d/2 

9M ^ d) = X 1/AoTi-l + d/a), d/2<a<d- (61) 
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4-2. Magnetization 

We first evaluate the large-time behaviour of the magnetization, m(t) := (sk=o(t))/N, 
which indeed obeys the expected scaling form [5]. Using equations ( 1T61) . ( 1561) and ( 1591) . 
we obtain the magnetization as 

-W = ^-^fl + -V 1/2 . (62) 



For t <C tM, the function m(t) ~ t e , and hence the initial-slip exponent is equal to 
9. While for t tM the asymptotic behaviour of the magnetization is m{t) ~ t e ~ 1 / 2 . 
Now the large-time behaviour of the order-parameter near the critical temperature is 
expected to relax as t~ /3 / 1/2 , where (3 and v are the equilibrium exponents, and the 
dynamic exponent z = a. For the long-range spherical model [12] (3 = 1/2 and 

r l/a, o <( x<rf/2 

1 l/(d-a), d/2 < a < d K J 

Thus in this model the slip exponent is not independent and is given by the relation 
6 = —(3/vz + 1/2. If we also introduce long-range correlations in the magnetized initial 
states, then 6 can depend on the classes of initial correlations, and hence can be an 
independent exponent. 

4-3. Auto- correlation and auto-response functions 

Since the expectation value of the spin is non-zero in the presence of initial 
magnetization, the Fourier transform of the connected correlation is given by 

C k (t, t w ) = C k (t, t w ) - N5 kj0 m(t)m(t w ). (64) 

Therefore, auto-correlation function, C(t,t w ) := C x= o(t,t w ) = J2k Ck(t, t w )/N, is given 
by 



C(t, t w 



mi)f[^\ + 2T c drfT—^-T g(r) .(65) 



y/g(t)g{t w ) l \ * J Jo 

Using the asymptotic form of g(t), as given in equation fl59l) . and f(t) ~ A t~ d / a ', and 
the fact that the above defined 9 < 1/2, equation fl65|) reduces to 

C(t,t w ) « A 2 d ^ tl- d ^F c (x,u), (66) 

where 

F c {x,u) = * lcX = dv{x + l-2v)- d/a v- 29 {l + uv) 1 (67) 

a/(1 + ux)(l + u) Jo 

x = t/t w > 1, and u = t w jtu- Hence the value of the exponent b = —1 -\- d/a. In the 
ageing regime, i>1, the above function takes the asymptotic form 

F c ( x , u) ~ — 2Tc ^ " ( + u—!—) , (68) 
K ' ^/(l + ux)(l + u) \l- 26 2 - 26 J' 1 ; 



The role of initial conditions in the ageing of the long-range spherical model 



13 



Conditions 


a = b 


Ac — 


X°° 


t w < t M , < a < d/2 


-1 + d/a 


d 


1/2 


t w < t M , d/2 < a < d 


-1 + d/a 


3d/2 - a 


1 - d/a 


t w >t M , < a < d/2 


-1 + d/a 


d + a/2 


2/3 


t w > tM, d/2 < a < d 


-1 + d/a 


(3d-a)/2 


d/(d + a) 



Table 3. Ageing exponents and asymptotic FDR in the Gaussian regime in presence 
of initial magnetization. 



and has ageing exponent Ac, as listed in the table®, which depends on whether t w <C tM 
or t w > t M . 

The auto-response function for the case with non-zero initial magnetization, from 
equations fl37j) and fl59l) . is obtained as 

R(t, t w ) = 2 d /°A t~ d ^F R (x, u), (69) 

where 



F R (x, u) = J ±±± (x - l)- d ^x e . (70) 
V 1 + ux 

Hence the value of the exponent a = —1 + d/a. The ageing exponent \r, which again 
depends on whether t w t M or t w ^> t M , is given in table®. Both the exponents, a 
and b, do not notice the time-scale tM and are unaffected by the initial magnetization. 
The expression for FDR, obtained from equations fl66l) and (169]) . reduces to 

X (tt\- T c F R (x,u) 

1 ' w) (1 - d/a + ud u - xd x ) F (x, uY (n) 
and results in the following asymptotic value of FDR, 

f (l-20)/(2-20), t u) «t M 
" \ 2(1- C)/(3 -2fl), t w ^>t M ' 1 ] 

which is explicitly given in the table®. 

In the case where t w -C tM, the initial magnetization is irrelevant and the ageing 
exponents are the same as those in the case of the uncorrelated initial state with no 
initial magnetization. As the waiting time t w increases and becomes larger than tM, the 
system will crossover to a new magnetized ageing class. 



5. Conclusion 



To summarize, the correlation and response functions for the non-conserved long-range 
spherical model are determined for two different sets of initial conditions: 

(i) long-range initial correlations with an arbitrary power a~o > — d, and vanishing 
magnetization; 

(ii) non-zero initial magnetization m , but no initial correlations beyond spherical 
fluctuations. 
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We evaluated the ageing exponents corresponding to various ageing classes. The ageing 
exponents can depend on the class of the initial state, though not on the details of it. 

When the long-range initial conditions are considered, the correlations and 
responses do show scaling behaviour for all cr > —d, in both the mean-field and the 
non-trivial cases of the long-range spherical model. 

In the mean-field case (0 < a < d/2), as the value of cr increases we distinguish 
three regimes: 

(i) regime II, where —d < o-q < a — d, 

(ii) regime V, where o~ > a — d and cr < — er, 

(iii) regime IV, where <To > o — d and a > —a. 

In both regimes II and V, the initial correlations are stronger than the equilibrium 
correlations, while in regime IV it is the other way around. We find that the exponents 
in regimes II and V can depend on do, whereas not in regime IV. In other words, we find 
that as we reduce the initial correlations there is a crossover from the initial-condition 
dependent ageing classes to that with no dependence. 

In the non-trivial case (d/2 < a < d), as the value of o~ increases we distinguish 
two regimes: 

(i) Regime I, where —d < cr < a — d, 
(i) Regime III, where o"o > o — d. 

We find that the exponents in regime I can depend on o"o, whereas this is not the case 
in regime III. Here too there is a crossover from the initial-condition dependent ageing 
class to that with no dependence, as <To goes from below to above a — d. In regime III, 
the initial correlations are weaker than the equilibrium correlations. While in regime I, 
depending on the value of a, it can be weaker, similar, or stronger, though the exponents 
depend of a . We lack a simple physical explanation of this dependence on the initial 
conditions in regime I. 

In the case of non- vanishing initial magnetization too, the correlations and responses 
show scaling behaviour. The initial magnetization m introduces a time scale tu-, and 
when the waiting time t w <C tu-, the ageing classes are same as those obtained with 
m = 0. As t w increases beyond tu-, there is a cross over to the magnetized ageing 
classes. 

Though we have evaluated only auto-correlation and auto-response functions, a 
similar analysis will suggest that the form of the two-point functions, in the presence of 
initial correlations, gets modified from that given in equation (CQ) to the following form, 

(0 4 (x,t)0 J (y,t w )> =C*F„ ( l -^jP-,f) +t~^F ij ( l -^jp-,f) ,(73) 

V t w l w J \ t w l w J 

where the second term on the right carries the memory of initial conditions, and defines 
new exponent a,j and scaling function Fij. The function asymptotically behaves as 
Fij(0,x) ~ x~ XiJ / z , and defines another exponent Ay. When Oy < a^-, then the 
fluctuations at time t within regions of length-scale much less than t 1 / 2 are no longer 
canonical equilibrium fluctuations. The above scaling form for the two-point functions 
is likely to hold beyond the model that we considered here. 
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